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Uncovering multiple Wolf-Rayet star-clusters and the ionized ISM 
in Mrk 178: the closest metal-poor Wolf-Rayet HII galaxy 
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ABSTRACT 

New integral field spectroscopy (IFS) has been obtained for the nearby metal-poor Wolf-Rayet 
(WR) galaxy Mrk 178 to examine the spatial correlation between its WR stars and the neigh- 
bouring ionized interstellar medium (ISM). The strength of the broad WR features and its low 
metallicity make Mrk 178 an intriguing object. We have detected the blue and red WR bumps 
in different locations across the field-of-view (~ 300 pc x 230 pc) in Mrk 178. The study 
of the WR content has been extended, for the first time, beyond its brightest star-forming 
knot uncovering new WR star-clusters. Using SMC/LMC-template WR stars we empirically 
estimate a minimum of ~ 20 WR stars within the region sampled. Maps of the spatial distri- 
bution of the emission-lines and of the physical-chemical properties of the ionized ISM have 
been created and analyzed. Here we refine the statistical methodology by Perez-Montero et 
al. (2011) to probe the presence of variations in the ISM properties. An error-weighted mean 
of 12+log(0/H)=7.72 + 0.01 is taken as the representative oxygen abundance for Mrk 178. A 
localized N and He enrichment, spatially correlated with WR stars, is suggested by this anal- 
ysis. Nebular Hen/i4686 emission is shown to be spatially extended reaching well beyond the 
location of the WR stars. This spatial offset between WRs and Hen emission can be explained 
based on the mechanical energy input into the ISM by the WR star winds, and does not rule 
out WR stars as the Hen ionisation source. We study systematic aperture effects on the detec- 
tion and measurement of the WR features, using SDSS spectra combined with the power of 
IFS. In this regard, the importance of targeting low metallicity nearby systems is discussed. 
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1 INTRODUCTION 

HII galaxies are local, gas-rich, systems with obsei^- 
able properties do minated by young and ma ssive star 
clusters (e.g., Melruck. Terlevich. & Eggletorj 1985c 



Campos-Aguilar. Moles. & Masegosal Il993l: iPopescu & Hopd 
20od iKehrig, Telles, & Cuisinieiiliodi : IWestera et al.ll2004l) . HII 
galaxies are also low metallicity objects [7.0 < 12+log(0/H) < 
8.3] in comparison with the solar metallicity [12+log(0/H)Q=8.69; 
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lAsplund et al.l ( 120091) 1 . and actually include t he most metal poor 
systems known in the local universe (e.g., iKehrig et al.l 120061 : 
IWesteraetal .ll2012l). for instance IZw 18 with 12+log(0/H) ^ 7.2 
(e.g.. lVilch ez & Igles ias-Pariunoll 19981 : IPapaderos & OstlirJlioTii 
and SBS 335-052W for which 12+log(0/H) w 7.0-7.2 (e.g., 
Ilzotov et al . 2009). Local metal-poor objects can be considered as 
template systems that help us to understand distant star-forming 
galaxies which cannot be studied to the same depth and angular 
resolution. 

Wolf-Rayet (WR) signatures (most commonly a broad fea- 
ture centered at ~ 4680 A or blue bump), indicating the pres- 
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ence of WR stars (the last observable stage during the evolution 
of massiv e stars before t he formation of neutron stars or black 
holes; e.g. lCrowthen '2007). have been found in the spectra o f some 
HII galaxies (e.g., JKunth & Sargent 1981; Legrand et aj 1 19971 : 
iGuseva. Izotov, & Thuanll2000l ; |Perez-Montero et al.l201ol) 7This is 
an important observational fact since single star, non-rotating stel- 
lar evolution models fail in reproducing the WR content in low 
metallicity environments (see 'Srinchmann. Kunth. & Durret 20081 
and references therein). The investigation of the WR content in 
galaxies is crucial to test stellar evolutionary models, specially at 
low metallicities where more data are needed to constrain these 
models. 

This work is part of our program to investigate HII galax- 
ies with WR features using integral field spectroscopy (IPS; 
iKehrig et alJl2008l ; |Perez-Montero et al]|201 ll) . IPS has many ben- 
efits in a study of this kind. Long-slit observations may fail in de- 
tecting WR features due to their faintness with respect to the stellar 
continuum emission and spatial distribution of WR stars across the 
galaxy. In particular in low metallicity objects, like HII galaxies, 
the dilution of WR features and the difficulty in spectroscopically 
identifying WR stars is even stronger owing to the steeper metallic- 
ity dependence of WR star wind s which lowers t he line luminosities 
of such stars (e.g. JCrowther & H adfield 2006). Kehrig et al.l ( l2008l) 
demonstrated for the first time the power of IPS in minimizing the 
WR bump dilution and finding WR stars in ex tragalactic systems 
where they were not dete cted before (see also ICairos et al.lCOld ; 
iGarcia-Benito et alj20Ioh . 

This technique also lowers the difficulty when studying the 
spatial correlation between WR stars and properties of the sur- 
rounding interstellar medium (ISM; e.g., Legrand et al. 1997; de 
Mello et al. 1998). Massive stars play a key role in the evolution of 
the ISM. The mechanical energy injected by evolving stars through 
wind s and supernova (SNe) explosions can re ach values ~ 10^^ erg 
(e.g., lMoIla. Garcia- Vargas, & Bressanll2009l) . The spatial distribu- 
tion of gas near the massive star clusters hosting WR stars may 
be affected by the disruption of the ISM, which also determines 
how newly synthesized metals are dispersed and mixed with the 
original gas from which the stars formed. Thus, a two-dimensional 
analysis of the ionized material in HII galaxies helps us to bet- 
ter understand the interplay between the massive stellar popula- 
tion and the ISM. Por instance whether WR stars are a significant 
contributor to abundance fluctuations on timescales of f ~ 10' yr 
and to the formation of hi gh-ionization lines (e.g. HeII/t4686) are 
stiU unsolved issues (e g., Roy & K unthll 19951 ; iKehrig et aJlUoill ; 
IShirazi & BrinchmannI 1201 2) that can be probed more precisely 
when applying IPS to nearby galaxies like Mrk 178, as we will 
show in this manuscript. The investigation on formation and there- 
abouts of gamma-ray bursts and Type Ib/c supernova progenitors, 
believed to be WR stars in metal-poor galaxies, may also bene- 
fit from the study presented here (e.g.. iWooslev & Bloomll2006t 
iModiaz et ani2008l) . 

The presence of a broadened 4686 A feature, attributable 
to WR stars, in the int egrated spectrum of Mrk 178 had 
alread y been noticed by iGonzalez-Riestra. Rego. & Zamoranol 
( Il988r). This system also appears in the WR galaxy catalogue 
by ISchaerer. Contini. & Pindaol 1 119991) and is contained in some 
studies of WR gal axies based on long-slit spect roscopy (e.g.. 



Table 1. General Properties of Mrk 178 



Parameter 



Mrk 178° 



Guseva. Izotov. & Thuan 2000; Buckalew , Kobulni ckv. & Dufoun 
2005l ; lBrinchmann. Kunth. & Durrell2008r) . In this work we present 



Other designation 

Morphological type* 

R.A. (J2000.0) 

DEC. (J2000.0) 

redshift 

D"(Mpc) 

Scale (pc///) 

(B-R)'' 

u'' (mag) 

g'' (mag) 

r^ (mag) 

f (mag) 

z' (mag) 

A,/^(mag) 

R2j^(arcmin) 

SFR(HQ')''(Moyr-') 

S5fR'(Mo yr-l Kpc-2) 



UGC 06541 

iI,M dwarf 

llh 33m 28.9s 

+49d 14' 14" 

0.0008 

3.9 

19 

0.89 ±0.10 

15.133 ±0.040 

14.228 ± 0.024 

14.139 ±0.023 

14.131 ±0.023 

14.153 ±0.021 

0.049 

0.62 

0.0085 

0.006 



" Mrk 178 belongs to the Ma rkarian l ists of galaxies w ith strong ultraviolet 
continuum (.Markari^l 19691) ; ' From lGil de Paz et al.l fcoOS): ^ Distance 
based on the ti p of the Red Giant Branch dis tance from lKarachentsev et alj 
)2O03h: '' FromlGil de Paz & Madord b005l) : •■ SDSS pho tometry from 
IShirazi & BrinchmannI j2012l) : f Galactic extinction from lSchlegel et al.l 
1199811 : ^ The optical radius R25 (= D25/2) in arcminutes, from the RC3 
catalog JD e Vaucouleurs et al. 1991); '' The Htv star formation rate from 
iLee et alj <2'o09l) : ' Star formation rate per unit area [= SFR/(;r X (R25)')]. 



the first two-dimensional spectroscopic study of the close-by HII 
galax y Mrk 178. At a distance of 3.9 Mpc jKarachentsev et alj 
l2003h and with a metallicity ~ I/IOZq, Mrk 178 is one of the most 



metal-poor nearby WR galaxies, which makes it an ideal object for 
our project. 

We focus on investigating the spatial distribution of nebular 
physical-chemical properties (e.g. electron temperature, gaseous 
metal abundances) and excitation sources for the gas in Mrk 178, 
and how these are related with the WR content derived here. 
We discuss the significance of the observed spatial variations 
of electron temperatures and chemical abundances found here, 
and the origin of the nebular Heii/14686 emission. Additionally, 
JBrinchmann. Kunth. & Durret! (|2008) found Mrk 178 to be an in- 
triguing object and significant outlier among their sam ple of WR 
galaxies from the Sloan Digital Sky Survey (SDSS; JYork et al.l 
2000). Based on the SDSS spectruin of Mrk 178, the authors mea- 
sured a remarkably high L(WR bump)/L(IIy6) ~ 0.05. This corre- 
sponds to a N(WR)/N(0) ~ 0.5 which is much higher than expected 
from any model at the low metallicity of Mrk 178. They suggested 
that this could be caused by the limited aperture size probed by the 
SDSS spectrum but were unable to test this possibility. In this work 
we test this hypothesis using our integral field unit (IPU) data. Ta- 
ble[T]compiles general properties of Mrk 178. In Pig. [T]we show a 
three-color broad band image from the SDSS of Mrk 178, and the 
HST/WPC3 image of Mrk 178 in the continuum filter P547M. 

The paper is organized as follows. In Sect. 2, we report obser- 
vations and data reduction. The WR content and nebular properties 
across the sampled regions of Mrk 178 are presented and discussed 
in Sect. 3. Pinally, Sect. 4 summarizes the main conclusions derived 
from this work. 
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Figure 1. Top: Three-band {z,r,i) SDSS color-composite image of Mrk 178 
overlaid with the observed field of view (FOV) of WHT/INTEGRAL (~ 
16 " X 12") represented by the white box. The Ha contours from Gil de 
Paz, Madore, & Pevunova (2003) are also shown. Bottom: HSTAVFC3 
image of Mrk 178 in the continuum filter F547M, coiTesponding to 
WHT/INTEGRAL FOV. The big (green) circle represents the approximate 
location of the SDSS fiber on Mrk 178. The small (blue) circle illustrates an 
individual spaxel (spatial element on the IFU) from the INTEGRAL fiber 
system. The three WR knots defined in section 3 are also marked. North is 
up and east is to the left. 



2 OBSERVATIONS AND DATA REDUCTION 

2.1 Observations 

IFS data of Mrk 178 wer e obtained with the INTEGRAL fiber sys- 
tem dArribas et al. 19980 in co mbination with the WYFFOS spec- 
trograph dBingham et al.ll994l) at the 4.2 m William Herschel Tele- 
cope (WHT), Roque de los Muchachos observatory. The data were 
taken with the 1024x1024 Tek6 CCD (Noise = 7 e; Gain = 2.1 
e/ADU). The observations were performed on 2001 January 18 un- 
der photometric conditions and atmospheric extinction in r' of ~ 
0.076 mag. 

We used the R600B grating which covers from ~ 37 10 to 6750 
A and provides a full width at half-maximum (FWHM) spectral res- 
olution of ~ 6A and a linear dispersion ~ 3A/pixel. The data were 
taken with the standard bundle 2, for which the field-of-view (FOV) 
is formed by 189 fibers (each ~ 0.9" in diameter) covering an area 



of 16" X 12" (~ 300 pc X 230 pc at the distance of 3.9 Mpc) on 
the sky (see Fig. [T}. One pointing of Mrk 178, encompassing its 
major star forming region towards the south, was observed during 
1.5 hours split into three exposures of 1800 seconds each. The sci- 
ence frames were taken at airmasses < 1.1 in order to avoid strong 
efi"ects due to differential atmospheric refraction. 

Additionally, calibration images (exposures of arc lamps and 
of continuum lamps) were obtained. Observations of the spec- 
trophotometric standard star Feige 34 were obtained during the ob- 
serving night to flux calibrate the data. 



2.2 Data Reduction 

The first steps of the data reduction were done through the P3d tool 
which was de veloped to deal with integral-field spectroscopy data 
dSandin et alj|201 0). After subtracting the bias level, the expected 
locations of the spectra were traced on a continuum-lamp expo- 
sure. The spectra were wavelength calibrated using the exposures 
of CuAr+CuNe arc lamps obtained immediately after the science 
exposures. We checked the accuracy of the wavelength calibration 
by measuring the central wavelength of the [OI]/15577 A sky line 
in all fibers and found a standard deviation of + 0.20 A. 

Fibers have different transmissions that may depend on the 
wavelength. The continuum-lamp exposures were used to deter- 
mine the differences in the fiber-to-fiber transmission and to ob- 
tain a normalized fiber-flat image, including the wavelength depen- 
dence. In order to homogenize the response of all the fibers, we 
divided our wavelength-calibrated science images by the normal- 
ized fiber-flat image. 

In the next step, the three exposures taken for the same point- 
ing were combined in order to remove cosmic rays, using the im- 
coMBiNE routine in IRAPU Flux calibration was performed using 
the IRAF tasks standard, sensfunc and calibrate. We co-added the 
spectra of the central fibers of the standard star exposure to create a 
ID spectrum that was used to obtain the sensitivity function. This 
sensitivity function was applied to the combined Mrk 178 science 
frame. 



3 RESULTS AND DISCUSSION 

3.1 Number and location of WR stars 

WR stars represent a late evolutionary phase of massive O stars 
whose initial mass (M,„;) on the main sequence was, at solar 
metallicity, > 25 Mq. At lower metallicities, this M„„ limit in- 
creases (e.g. Crowther 2007). According to de Mello et al. (1998), 
the M,„, is > 90 M© for Z ~ 1/30 Zq, based on non-rotating 
models for single stars. Nevertheless, a significant number of 
WR stars detected in very metal-poor objects appears contrary 
to the expectations of such models (e.g. Izotov et al. 1997; Pa- 
paderos et al. 2006). Including rotation in the massive stellar evo- 
lution lowers the minimum mass required for a star to enter the 
WR phase, increasing their population in metal-poor environments 
dMevnet & Maeder 2005). WR stars are also extremely hot and lu- 
minous with temperatures that sometimes exceed 50,000 K and 
some with luminosities approaching 10* L© (iCrowther. .2007.) . The 



' IRAF is distributed by the National Optical Astronomical Observato- 
ries, which are operated by the Association of Universities for Research 
in Astronomy, Inc., under cooperative agreement with the National Science 
Foundation 
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Figure 3. From top to bottom: dereddened spectra of knots A, B, and C (black lines), corrected for a distance of 3.9 Mpc. Generic blue and red WR features 
[WN templates (red lines); WC/WO templates (green lines); composite WN+WO/WC templates (blue lines)] are shown for comparison. Templates of SMC 
WR stars are used where possible (WN,WO). We adopt WC4 LMC template since no SMC counterparts are available. The spectra of the WR knots and WR 
star templates are continuum subtracted. 
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Figure 2. Portion (~ 4300 - 6000 A) of the integrated spectrum, in units of 
erg s"' cm"^ A"' , for the three knots in which the WR features are detected 
(see Figs.|5][2]and[8). The spectral range for both blue and red WR bumps 
are marked. 



spectra of WR stars show strong, broad emission lines formed in 
their fast, dense winds. Two main classes of WR stars are defined 
based on emission line ratios: nit rogen-type (WN) and carbon- 
type (WO and WC) stars (e .g., 'Smit h, Shara, & Moffall 1 19961 : 
ICrowther, De Marco, & Barlow 1998). The so-called WN types are 
nitrogen-rich heliu m stars that sh ow the results of the CNO hydro- 
gen burning cycle ( iSmithlll973h . Besides the classic He-burning 
WN stars, we also find in the liter ature examples of H-rich WR 
stars usually very luminous(WNh; ISmith. Shara. & Moffka 1 19961 : 



iFoellmi. Moff^at, & Guerrerol 120031: ISchnurr et ahlbOOSl) w hich are 
believed to be core-H burning (e.g. ICrowther et al.lll995h . Neither 
H nor N are detected in the spectra of carbon-type stars which show 
predominantly the products of the triple-o- and other helium burn- 
ing reactions (see Crowther 2007). 

Integrated spectra of star-forming galaxies sometimes reveal 
WR spectral features, called WR bumps, which allow WR stars to 
be identified in both Local Group and more distant systems. Ex- 
haustive searches for extragalactic WR stars have been the sub- 
ject of several works (e.g., Buckalew, Kobulnicky, & Dufour 20051 ; 
ICrowther et"Zll2007l; iNeugent & Massevll201lh . The blue md red 
WR bumps observed in galaxy spectra are centered at ~ 4680 and 
5808 A, respectively. While WN stars are mainly responsible for 
the formation of the blue bump, the carbon-type WR stars con- 
tribute the most for the red bump emission. 

WR features in the s pectrum of Mrk 178 have been re- 
porte d in previous work fOonzalez-Riestra, Rego, & Zamorarid 
i l98^ ; ISchaerer, Contini, & Pindao 1999; Guseva, Izotov, & Thuan| 
1 20001) . From our IFU data we detect WR features in three differ- 
ent regions in the galaxy which we defined as WR knots (the lo- 
cation of the WR knots are marked on all maps of line intensity, 
line ratio and ISM properties; see Figs. |5]|7]and[8}. Knots A and 
B are composed by three spaxel^ while knot C is represented by 
a single spaxel. The detection of WR signatures at the locations of 
knots A and C is reported for the first time in this work. We created 



^ Individual elements of IFUs are often called 'spatial pixels' (commonly 
shortened to 'spaxel'); the term is used to differentiate between a spatial 
element on the IFU and a pixel on the detector 
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one-dimensional spectra for knots A and B by adding the emission 
from the coiTesponding spaxels. Portions of the resulting spectra 
are shown in Fig.|2]from which we can see the WR signatures. The 
spectrum of knot A only displays the blue WR bump while we see 
both blue and red WR bumps in knots B and C. No WR features 
were detected in any parts other than Knots A, B and C across the 
sampled regions. We shall now discuss the WR content in Mrk 178 
within the INTEGRAL FOV. 

According to Starburst 99 synthesis evolutionary models (Lei- 
therer et al. 1999) the expected number of WR stars for a star 
cluster at an age of ~ 4 Myr (the age at which models predict 
the WR content to be maximum), and with the Ho- luminosity 
(~ 10"'^ erg/s) and metallicity similar to those measured for the 
three WR knots in Mrkl78 is < 1. State of the art models of stel- 
lar evolution that include the effects of rotation or binary evolu- 
tion on the main sequence and beyo nd might be able to repro- 
duce the observed WR features (e.g.. lEldridge & Stanwavll2009l 
and references therein). However, since full sets of tracks are not 
yet available and star forming regions do not fully sample the up- 
per IMF, a more empirical approach has been taken. In this re- 
gard, we will use template fitting of individual WR stars to es- 
timate the WR star content in this work. This methodology has 
been successfully applied in synthesizing the W R bumps in nearby 
spiral and irregular galaxies (e.g.,ISid oli. Smith . & Crowthen2006l : 
iHadfield & Crowtheilliooel; IMoU et al..,2007,) including one exam- 
ple for a metal-poor [I2+ log(0/H) < 8.0] HII galaxy (IZwl8; 
ICrowther & Hadfieldll2006t) . 

Here, to assess the number of WR stars in each knot, based 
upon c urrently available templ ates, we use SMC [12+log(0/H) 
K 8.0; iRussell & Dopit j ( 119901 1 template WR stars [early-type 
WN (W N2-4), mi d-type WN (WN5-6) fromlCrowther & Hadfieldl 
( l2006l) : WO from lKingsburgh. Barlow. & StorevI l ll995l) 1 because 
of the low metallicity of Mrk 178 [12+log(0/H) » 7.72; see sec- 
tion [TZS). LMC templates are used solely when SMC templates 
are unavailable (e.g. WC subtypes). These template WR stars were 
adjusted for the redshift of Mrk 178 (250 km s"'). 

Figure |3] shows the comparison between the observed WR 
bumps for our 3 knots and the WR template spectra [predic- 
tions for WN stars (red), WC/WO stars (green) and composite 
WN+WC/WO stars (blue)]. In view of the limited templates avail- 
able and variations in line strength and line width for individual WR 
stars we wish to emphasise that the results presented here are not 
unique, but that we attempt to match the morphology (line strength 
and width) using the fewest number of subtypes. 

In Knot A, a broad blue (A^ -4686; FWHM~I5 A) emission 
feature is observed, together with possible nebular Heii/i4686 com- 
ponent. The FWHM (~5 A) of this narrow Hen component, compa- 
rable to that of the other nebular emission lines (e.g., [O iii]/15007), 
and its spatial extent are evidence of its nebular nature. Applying 
the SMC WN calibration from lCrowther & Hadfieldl ( |2006|) to the 
blue bump would require ~ five equivalent WN5-6 stars. 

For knot B, we estimate that six WN5-6 and eight WO stars 
[CIV/i5808 is very broad (FWHM- 100 A) requiring a WO instead 
of WC template] are needed to reproduce the WR features. The 
blue WR spectral feature in knot C is fairly well reproduced by a 
population of 1 SMC-like WN5-6 star and 1 LMC-like WC4 star 
(LMC template WC4 stars were used since no WC4 have been ob- 
served in the SMC). However the broad red feature requires a WO 
star to match the observed morphology, while the blue comparison 
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Figure 4. This figure illustrates the main output from spectral synthesis for 
a single spaxel (~ 1" from Knot B) in the region of Mrk 178 sampled in 
this work. The blue and red lines (top panel) con'espond to the observed 
spectrum and the best-fitting synthetic SED for the stellar continuum, cal- 
culated with the population synthesis code STARLIGHT, respectively. By 
subtracting the latter from the former we obtain the pure nebular emission 
line spectmm (green in the bottom panel). The y-axis is in 10"'^ erg s"' 
cm"- A" ' . The smaller panel shows the age distribution of the stellar popu- 
lations in the best-fitting synthetic SED as a function of their luminosity (at 
4200 A) in percent. 



is much worse. Alternatively, we checked whether a WN/C stajj 
might enable a match to both blue and red features simultaneously, 
but a satisfactory solution could not be achieved since both fea- 
tures are much too broad with respect to a WN/C star. While the 
small number of appropriate templates from the SMC means that 
we cannot achieve a perfect fit for knot C, we find that a population 
of ~ 1-2 WR stars is likely. 

Of course, one needs to bear in mind that the numbers of WR 
stars implied from these comparisons are only as good as the tem- 
plates used. For instance, there are only two carbon sequence WR 
stars at metallicities below the LMC in the Local Group: the WO 
stars DRl i n IC 1613 ( Kingsburgh & Barlow 1995) and Sand 1 in 
the SMC ( Kings burgh. Barlow, & Storevll995l) . In view of this, we 
have also examined the DRl template, but the width of the red 
bumps observed in knots B and C is in much better agreement with 
the SMC template. In addition, the use of LMC templates for WC 
stars (in the absence of SMC counterparts) is clearly not ideal. Still, 
the morphology of WC stars in the metal-poor dwarf IC 10 is rem- 
iniscent of the LMC. Indeed, the average C iv /1/1580I-I2 line lu- 
minosity of IC 10 WC stars is consistent with that of the LMC WC 
population, in spite of large uncertainties due to high foreground 
interstellar extinction (Crowther et al. 2003). 

Early-type WN templates would indicate a larger WR popula- 
tion (factor of ~ 10 times higher than for mid-type WN stars), al- 
though we prefer mid- type WN stars throughout to early- WN stars 
because they provide a better match to the Heii/14686 line width. 
Of course, individual WN stars at a given subtype span a range of 
line widths , but the range is modest for mid- and late-type stars (see 
Fig.2 from lCrowther & Hadfiel32006i) . 

In addition to the observed line width constraint, mid- 
WN stars tend to be especially common in giant HII regions 
across a range of metallicities, including our Milky Way (e.g., 

^ Our template, Brey 29 (BAT99-3 6), is the LMC composite WN4bAVCE 
JFoellmi. Moff'at. & Guerrerol2003l) 
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NGC 3603; Drissen et al 



Moffat, Seggewiss & Shars Im 



T995h . 
1981; 



the LMC0(e.g. 30 Doradus; 
and the SMC (e.g. NGC 346; 



Conti. Garmanv. & MassevI 1 1989 1. Therefore, we might expect 



these to dominate over early WN stars in knots A and B for which 
the spatial extent is ~ 4 arcsec (~ 80 pc at a distance of ~ 4 Mpc), 
the size of typical HII regions. The size of Knot C (radius ~ 8 
pc) is of the or der of some ejecta nebula e around isolated/binary 
WR stars (e.g., Esteban & vi lchezl 1 19921; IStock & Barlow' 2010; 
iFemandez-Martm et al. 2012) which supports our estimation of 
very few WR stars (~l-2) there. It is likely that the number of 
WR stars estimated here represents a lower limit owing to the use 
of SM C templates and reduced lin e luminosities at lower metal- 
licities JCrowther & Hadfieldll2006h . Unfortunately, a more quan- 
titative assessment would require more individual metal-poor WR 
stars to be included in low metallicity templates (ICIO is the only 
suitable candidate at the present time). 



3.2 A 2D view of the warm ISM 

3.2.1 Line fitting and map creation 

To derive line ratios properly, we first subtracted from the ob- 
served spectra the spectral energy distribution (SED) of the un- 
derlying stellar population found by the spectral s ynthesis code 
STARLIGHTJl icid Femandes et alj l20o4 120051; iMateus et al.l 
12006) . STARLIGHT fits an observed continuum SED using a com- 
bination of the synthesis spectra of different single stellar pop- 
ulations (SSPs) using a x^ minimiza tion procedure. W e chos e 
for our analysis the SSP spectra from Bruzual & CharlotI (120031) . 
based on the STELIB library of iLe Borgne et al.l ( |2003|k Padova 
1994 evolutionary tracks and a lchabrieij (l2003n initial mass func- 
tion between 0.1 and 100 Mq. We used the library metallicities 
(Z=0.00 1,0.004,0.008) closest to the average oxygen abundance 
measured in the gas [124-log(0/H) » 7.72; see Section [3.2.31 and 
a set of 38 ages from 1 Myr up to 10 Gyr The reddening law 
from ICardelli et al.l ( Il989l) with R^ = 3.1 was applied. Bad pix- 
els and emission lines were excluded from the final fits. According 
to the model SSPs, the correction factor for stellar absorption in the 
EW(H;3) is ~ 10% on average. Our nebular lines have EW(Hy8) > 
10 A so the effect of the underlying stellar population in Balmer 
lines involved in the reddening derivation is negligible. The best 
spectral synthesis fits also show that young stellar populations (age 
< 10^ yrs) dominate the light fraction in most spaxels. An illustra- 
tive example of the observed, modeled, and emission-line spectrum 
for a single spaxel in the FOV of Mrk 178 is displayed in Fig.|4] 

After subtracting the underlying stellar population for each 
spaxel spectrum, the emission-line fluxes are measured using the 
IRAF task splot. The intensity of each emission line was derived 
by integrating between two points given by the position of a lo- 
cal continuum placed by eye. The line-flux errors were calculated 
using the following expression: 



crime = 0-„„„N^'- I 1 -I- 



EW\ 
NAA 



1/2 






(1) 



where CTcont is the standard deviation of the continuum near the 
emission line, A' is the width of the region used to measure the 

'* Crowther & Hadfield (2006) included 15 LMC WN5-6 stars in their mid- 
WN template, of which 7 lie within 30 Doradus (R144-147,R136al,a2,a3). 
^ The STARLIGHT project is supported by the Brazilian agencies CNPq, 
CAPES and FAFESF and by the FranceBrazil CAFES/COFECUB pro- 
gramme. 



line in pixels, AA is the spectral dispersion in A/pixel, and EW rep- 
resen ts the equivalent width of the line JGonzalez-Delgado et al.l 
1 19941) . This expression takes into account the error in the contin- 
uum and the photon count statistics of the emission lineQ The rela- 
tive errors in the line intensities are ~ 5-8% for the bright lines (e.g. 
[Oii]/13727; [Oiii]/15007; Ha). Typical uncertainties for the weak- 
est lines ([Oiii]/14363; Heiii4686; [Nii]i6584) are about ~ 15 % 
and may reach up to ~ 20-30% in some fibers. The final errors do 
not account for flat-field and instrumental response uncertainties 
which are not significant (< 5%). 

Using our own IDL scripts we combine the flux intensities 
with the sky position of the fibers to create the maps of emis- 
sion lines and line-ratios presented in this paper. Figure |5]displays 
flux maps of some of the relevant emission lines and line ratios. 
As a guide to the reader, the spaxels where we detect WR sig- 
natures are indicated in all maps. The Ho- and [On]/i3727 maps 
show a larger area than the Hen/i4686 map because those lines 
are among the brightest optical emission lines in our data sets. 
The maps of [Oii]/13727/HyS and [Oiii]i5007/H8 present opposite 
trends; [Oii]/l3727/IIy8 peaks in the outer zones of the FOV while 
the highest values of [Om]/i5007/IIj8 are observed close to the WR 
knot thereabouts where the ionizing flux should be harder (see sec- 
tion 13.2.41 for discusson on excitation sources). The map of the 
[Oii]/13727/[Oiii]/i5007 ratio which is an indicator of the ioniza- 
tion parameter, is also shown. In addition, Fig|6] displays the Ho- 
map along with representative emission-line spectra from spaxels 
in different positions across our FOV. 



3.2.2 Physical conditions and chemical abundances in the ISM 

The reddening coefficient corresponding to each fiber spectrum, 
c(Hy3), was computed from the ratio of t he measu red-to-theoretical 
Hc/Hy? assuining the reddening law of ICardelli et al. (1989), and 
standard HII region characteristics (T^ = lO'* K; n^ = I00_cn£) 
which give an intrisic value of Ha/HyS = 2.86 (.Storev & Hummen 
1 19951) . 

Using the de-reddened line-ratios, we derived the physical 
properties and ionic abundances of the ionized gas for Mrk 178 
following the 5-level atom FIVEL program ( IShaw & DufouilI994l) 
available in the tasks TEMDEN and IONIC of the STSDAS pack- 
age. We calculated the final errors in the derived quantities by error 
propagation and taking into account eiTors in flux measurements. 
Systematic errors are not included. 

We obtained the electron densities, n^, from the 
[Sii]/16717//i6731 line ratio. The derived estimates and upper 
limits for n^ place all of spaxel spectra in the low-density regime 
(n, <300 cm^). 

We derived the T^ values of [Om] using the 
[Oiii]^4363/[Oin]i4959,5007 line ratio. r,[Oii] was calculated 
from the relation between [On] and [ Om] electron temperatures 
given bv lPerez-Montero & Diazl ( 120038 ) assuming an electron den- 
sity of 100 cm"^. We measured the faint auroral line [Oiii]/i4363 
for 60 spaxels with S/N > 3. These spaxels cover a projected 
area of nearly 42 arcsec' equivalent to ~ 0.02 kpc^, including the 
central starburst region. 

The oxygen ionic abundance ratios, O^/H^ and 0^*/H^, were 
derived from the [Oii]/13727 and [Om]AA 4959, 5007 lines, respec- 
tively using the corresponding electron temperatures. The total 



We caution that for faint stars near the detection limit CCD readout noise 
can also be important. 
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Figure 5. Intensity and line ratio maps of Mrk 178: Ua, narrow Heii/i4686, [Oii]^3727. [Oii]/i3727/H/3, [Oiii]/i5007/H/3 and [Oii]^3727/[Oiii]/i50()7. Dotted 
circles show the spaxels where no measurement is available. All maps ai'e presented in logaiithmic scale. The three WR knots (A, B and C) are labeled on the 
Ha flux map, and the spaxels where we detect WR features are marked with green crosses on all maps. For each map the vertical axis of the right-hand side 
color bars represents the values of the corresponding line (-ratio) while the horizontal axis (cr) displays the uncertainties associated with each quantity, which 
are coded using an ff-channel (transparency of the spaxel color). According to the ff-channel coding, the more intense the color the lower the uncertainty. North 
is up and east to the left. 



oxygen abundance is assumed to be: 0/H = 0^/H^ + O'^/H^ . With 
regards to the oxygen ionization correction factor (ICF), a small 
fraction of O/H is expected to be in the form of the O'* ion in high 
excitation HII regions when Heii/i4686 emission line is detected. 
In our case, we find the correction for the unseen higher ionization 
stages of oxygen to be negligible; in the spaxels where we measure 
Hen/14686, the ICF(O) correction would increase the total O/H by 
less than ~ 0.02 dex, which is considerably lower than the typical 
error associated with our O/H measurements (see top-left panel of 



Fig.|7j. N^ abundances were derived using the [Nii] emission line at 
6584 A and assumming rj[Nii] ~ r(,[Oii]; the N/0 abundance ratio 
was derived under the assumption that N/O = N^/O^. Finally, b ased 
on the Hei/16678 and equations from ;01ive & SkillmanI ( l2004h . we 
calculate He^/H^. These equations were used for the appropriate 
values of n^, temperature, and account for collisional deexcitation. 
We ignore Hei/15876 since it is affected by Galactic Nal absorption 
and Hei/14471 due to too low S/N. 

We also derived O/H values based on the N2 parame- 
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Figure 6. Representative spectra in units of 10 '^ erg s ' cm ^ A ', con'esponding to the nebular emission of three single spaxels across the FOV. 



ter [=log([Nn] ^6584/Hq')] by IPerez-Montero & Continil ( l2009l) 
to have an estimation of the metallicity over a larger number 
of spaxels. The N2 parameter is almost independent of redden- 
ing effects, and the relationship between N2 and O/H is single- 
valued. Moreover, to derive N/0 in a larger spatial extent, we 
made use of the N202 parameter [= log([Nn].16584/[On]^3727)] 
which has a linear relation with N/O at all ranges of metallicity 
jPerez-Montero & Continill2009h . 

Maps and histograms of the derived O/H and N/O are dis- 
played in Figs.|7]and[8l respectively. 

For the three WR knots A, B and C, the physical-chemical 
properties are derived following the methodology described above 
and are presented in Table |2] 



3.2.3 Spatial variations of the ISM properties 

Here, we refine the method described in ' Perez-Montero et al.l 
( 1201 If) to study spatial variations of ISM properties based on statis- 
tical criteria. We assume that a certain physical-chemical property 
is homogeneous across our INTEGRAL FOV if two conditions are 
satisfied: for the corresponding dataset (i) the null hypothesis (i.e. 
the data co me from a nor mally distributed population) of the Lil- 
liefors test ( lLillieforslll967.) cannot be rejected at the 10% signifi- 
cance level, and (ii) the observed variations of the data distribution 
around the single mean value can be explained by random errors; 
i.e. the coiTesponding Gaussian sigma (cr amission) should be of the 
order of the typical uncertainty of the considered property; we take 
as typical uncertainty the square root of the weighted sample vari- 
ance {cr„gigh,gj). In Table[3]we show the results from our statistical 
analysis for relevant ISM properties. 

For the electron temperature, despite the fact that the distri- 
bution of its measured values can be represented by a Gaussian 
fit (based on the Lilliefors test), we find that a-omissian > CTtveigiued 
(see Table [3}. This suggests that as a first approximation ran- 
dom variables alone should not explain the T^ distribution. The 



source of these variations is hard to assess. A possible explana- 
tion could be some error source that we do not take into account 
(e.g. variable seeing). However if this was the explanation, this 
source of error should affect all emission lines which does not 
seem to be the case. Our statistical study indicates that the scat- 
ter in T,, can be larger than the scatter in O/H within Mrk 178. 
IVilchez & Iglesias-Paramd 1 119981) find similar results for IZwl8. 
The opposite is also found in the literature; i.e. chemically homo- 
geneous gas is found in star-forming regions with no apparent T, 
variation (e.g.. lGonzalez-Delgado et al.lll99q) . 

From Table [3] we find that the measured values of O/H, both 
from direct-method and N2 calibrator, are fitted by a normal dis- 
tribution according to the Lilliefors test; additionally a-oaussian ~ 
(^weighted in both methods. We find the same results when consid- 
ering only the O/H values from the spaxels where we measure the 
nebular Hen/14686 (gray bars in the histograms of Fig. |7). Thus, 
according to our statistical analysis, the oxygen abundance in the 
ISM of Mrk 178 is homogeneous over spatial scales of hundreds of 
parsecs. However, we must bear in mind that this statistical method- 
ology ignores any spatial information and cannot be used to discard 
the presence of small-scale localized (~ 20 pc, our resolution ele- 
ment size) chemical variations, as we will discuss further in this 
section. 

Close inspection of the two histograms in Fig. |7] indicates 
that the Gaussian fit for the distribution of O/Hni is slightly 
broader than that for the OfHre histogram ((To„„j„„„,v2=0. 16; 
o"Goi/MiaH,rf=0.12), and the Gaussian peak is somewhat shifted to- 
wards higher values. This probably reflects the dependence of N2 
on the ionization parameter. Still, the mean values of the corre- 
sponding 0/Hjv2 and O/H7-J distributions are consistent within the 
errors (see Table|3)- 

In the case of the N/O ratio, the width and peak of the his- 
tograms for both (N/0)Te and (N/0)a>202 values are similar (see Ta- 
ble[3]and Fig.[8}. The distribution of N/O values can be represented 
by a Gaussian fit in agreement with the Lilliefors test, indepen- 
dently of the method used to estimate the N/O ratio. Also, a^Gaussimi 
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Figure 7. Top panels: maps of oxygen abundance (12+log O/H) derived from Tg[Om] (left) and the N2 parameter (right); right-side color bars and green 
crosses are as indicated in Fig. [5] The conesponding histograms are in the bottom row. The gray bars indicate the spaxels where we measure the nebular 
Heii/i4686, blue bars represent the spaxels showing only the blue WR bump (i.e. the spaxels within knot A; see Figs.|2]and|5), red bars indicate the spaxels 
within Knot B (see Figs.[2]and[5l, all showing both the blue and red WR bumps, and the black bar represents the individual spaxel in the NW with both the 
blue and red WR bump detection (i.e. Knot C in Figs.[2]and[5). The solid and dashed curves are the Gaussian fits to all spaxels and only the spaxels where 
we detect nebular Hen/i4686 line, respectively. The uncertainties of (0/H)a;2 are obtained by propagating the errors of [Nn] and Ho- lines (involved in the N2 
parameter) in the derivation of O/H, and do not consider the typical uncertainty associated with the N2 parameter (~ 0.3 dex). 



is of the order of cr„,p,j,,,„^ as is the case with O/H. Taking into ac- 
count only the spaxels with nebular Heii/14686 detection (gray bars 
in the histograms of Fig. [8j yields the same results. So, overall, 
across the whole FOV (~ 300 pc x 230 pc) of Mrk 178 we find no 
statistically significant variations in N/O. 

Since our statistical analysis does not consider any particu- 
lar spatial information, as mentioned above, to probe the existence 
of small-scale (~ 20 pc) localized chemical pollution, we have 
inspected the chemical abundances for the WR knots as derived 
in section 3.2.2 (see Table 2). From Tableland the N/O distri- 
butions in Fig. [8] we note that WR knot C appears to present a 
higher N/O. Since N enrichment fro m WR star winds is expected 
to come with He enhancement (e.g.,|Pagel. Terlevich, &Melnickl 
1 19861 : lEsteban&Vilchezll 19921 : iKobulnickv & Skillmarj |l996'), we 
analyse the He*/H^ as well. From Fig.|9l the histogram of He^/H^, 
and Table [2] one sees that WR knot C presents a considerably high 
He^/H+ (~ 0.18) in comparison to the rest of the spaxels. Thus, 
the existence of a localized N and Hqj pollution by WR stars 

^ Considering the high ionization degree of the WR knots (O^^/O is > 0.70 



in knot C cannot be ruled out. Since iMaeden ( Il992h , it is well 
known that WR stars are significant contributors of He, N and C. 
At these small-scales, chemical enrich ment has been previously 
reported for Galactic WR nebulae (e.g.. lEsteban & Vilchezl ll992l : 
Femandez-Martin et al. 2012) and for irregular dwarf galaxies, for 
insta nce IC 10 ( iLopez-Sanchez et al. 2011* ) and NGC 5253 (e .g., 
iKobulnickv et aU1997, : ,Monreal-Ibero. Walsh. & VilchezI 20121) . A 
statistical detection of N en hancement in galaxies showing WR fea - 
tures has been reported bv iBrinchmarm. Kunth. & Dur ret (20081). 
No evidence of chemical enrichment in the gas associated with the 
WR knots A and B is found. Nevertheless, one should bear in mind 
the detection of any chemical pollution on scales of ~ few parsecs 
will escape this work because of our resolution element size (~ 20 
pc at the distance of Mrk 178). 

Summarizing, chemically homogeneous ionized gas on 
spatial scales of > 100 p c appears to be a rul e in HII galaxies (e.g., 
iLee & Skilliiiaril |2004| : iKehrig et alj l2008l : iPerez-Montero et al.l 



for the three WR knots), their corr esponding ICF(He) is expected to be ~ 1 
ilzotov. Thuan. & Stasinskal2007h 
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Figure 8. Top panels: maps of log(N/0) derived from T^ [Oiii] (left) and the N202 pai'ameter (right); right-side color bars and green crosses are as indicated in 
Fig. \5\ The corresponding histograms are in the bottom row; the colouring of the bars and curves are the same as in Figure [7] The uncertainties of (N/0)jv202 
are obtained by propagating the errors of [Nil] and [On] lines (involved in the N202 parameter) in the derivation of N/0, and do not consider the typical 
uncertainty associated with the N202 parameter (~ 0.25 dex). 



20091; ICairos et all l2009allbl; |Perez-]V[ontero et al.l 1201 iL 
Garcia-Benito & Perez-tvlonterd l2012h . One may speculate 



that the ionized ISM in IVIrk 178 and other HII galaxies appears 
presently well homogenized because the freshly synthesized ele- 
ments were rapidly dispersed, on timescales < 10' yr, and mixed 
very quickly in the ISM. This hypothesis requires that ejected 
metals are both homogeneously dispersed and mixed with the 
ambient ISM on time scales comparable to that of the HII regions 
which seems to be unlikely b ased on hydrodynami c al and physical 
proce s ses in the gas (see iRov & KunthI 119951 : iTenorio-Tagld 
1 1994 iKobulnickv & SkillmanI Il997l) . Alternatively, the scenario 
proposed by IXenorio-Tagl e f 1996^ in which the ejecta from stellar 
winds and supemovae could undergo a long (~ 100 Myr) cycle in a 
hot phase (~ 10' K) before mixing with the surrounding ISM, has 
been discussed in the lite rature (e.g. IKobulnickv & SkillmanI 19971 : 
Ivan Zee & Havnej 120060 as a possible explanation for the homo- 
geneous chemical appearence of HII galaxies on spatial scales of 
> 100 pc. However the fate of the metals released by massive stars 
in HII regions is still an open question, and the proce sses of metal 
dispersal and mixing are very difficult to model (see iRecchi et al.l 
l2003l : lRecchi & Hensleill2013h . Thus, to better understand how the 
ejected metals cool and mix with the ISM, further investigation 



on the metal content of the differe nt ISM phases is needed but is 
beyond the scope of this work (e.g. lLebouteiller et al.l2009l) . 

In this work we find that the representative metallicity of 
Mrk 178 is 12H-log(0/H) = 7.72 ± 0.01 (~ 1/10 of the solar metal- 
licity) which represents the derived error-weighted mean value of 
0/H and its corresponding statistical error from all the individual 
spaxel measurements obtained from the direct-method (for this 0/H 
distribution, the square root of the weighted sample variance asso- 
ciated is 0.13 dex; see Table[3]l. 



3.2.4 Excitation sources and nebular HeuA4686 emission 

Although photoionization by massive stars is the most likely 
main excitation source in HII regions, some questions still re- 
main regarding the origin of high-ionization nebular lines as 
for instance Heii/i4686. This line has been observed in sev- 
eral WR gala xies, especially in the low-metal licity ones, like 
Mrk 178 (e.g., ISchaerer. Contini. & Pindaol [19991) . A possible ex- 
planation for this is that in such metal-poor environments, WR 
stars possess sufficiently weak winds that are optically thin in 
the He^ continuum (e.g.. lSchmutz. Leitherer, & Gruenwaldlll992l : 
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Table 2. De-reddened emission line-fluxes relative to 1000-I(H/3) and physical properties for the WR knots. The properties of the WR bumps are quoted. 



Wavelength 



3727 [O II] 
4363 [O III] 
4686 He II 
4861 H/3 
5007 [O III] 
6563 Ha 
6584 [N II] 
6678 Hel 
6717 [S 11] 
6731 [Sll] 



c(H/?) 

-EW(H/3) (A) 
logF(H8(erg/s) 



«e([S ll])(cm-J) 
r,([0 111]) 
12+log {0+/H+) 
n+log (0++///+) 
12+log {N+/H+) 
12+log (0/H)r. 
log (N/0)r, 
He+/H+ (i6678) 
He2+/H+ (i4686)° 



Knot A 



KnotB 



Knot C 



1320 ±60 


1670 ± 110 


955 ± 60 


135 ± 14 


112 ± 17 


147 ± 12 


40 ±4 


— 


— 


1000 ± 20 


1000 ± 90 


1000 ± 50 


5290 ± 160 


5213 ± 360 


5330 ± 180 


2810 ±60 


2810 ±240 


2770 ± 110 


51 ± 14 


93 ±30 


77 ± 16 


37 ±6 


— 


63 ±16 


112±5 


132 ± 11 


129 ± 15 


68 ±5 


89 ±9 


102 ± 16 



0.16 ±0.03 
80 ±2 
-14.00 



0.26 



0.10 
3± 1 
14.25 



0.14 ±0.05 
33 ±5 
-14.75 



<10 
16186 ±675 
7.10 ±0.08 
7.65 ± 0.05 
5.65 ±0.14 
7.76 ± 0.06 
-1.45 ±0.16 
0.11 ±0.02 
0.004 ± 0.001 



<10 
14600 ± 750 
7.29 ±0.12 
7.76 ± 0.07 
5.97 ±0.17 
7.89 ± 0.08 
-1.32 ±0.20 



210 
17079 ± 530 
6.93 ± 0.08 
7.60 ± 0.03 
5.81 ±0.10 
7.69 ± 0.09 
-1.12 ±0.13 
0.18 ±0.05 



log L(WR blue bump)* (erg/s) 


36.60 ± 0.07 


37.52 ± 0.06 


37.22 ± 0.07 


-EW (WR blue bump) (A) 


7±2 


8±2 


50 ± 12 


log L(WR red bump)* (erg/s) 


- 


37.88 ± 0.05 


36.87 ± 0.06 


-EW (WR red bump) (A) 


- 


20 ±6 


55 ± 15 



(a) estimated from Heii/i4686 and TeCLOiii]) llzotov. Thuan. & Lipovetskvil 19941) . (b) extinction-corrected luminosity of the WR bumps derived at our 
assumed distance of 3.9 Mpc; to measure the WR features we fitted the stellar continuum using a polynomical function and we subtracted it from the WR 
emission. Nebular Heii/14686 line was removed before deriving the luminosity and EW of WR blue bump for knot A. 
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Figure 9. Histogram of the He*/H^ ionic abundance ratio; the colouring of 
the bars and curves are the same as in Figure |7] 



ICrowther & Hadfieldl 120061) , therefore allowing the He^-ionizing 
photons to escape from the stellar atmospheres. 

However, WR features are not always seen when nebular Hen 
is observed. The nebular emission of Hen/14686 detected in this 
work is shown to be spatially extended, and some Hen-emitting 



regions are not coincident with the location of the WR bumps by 
tens of parsecs (see the map of Hen/14686 in Fig. O. Similar r e- 
sults have been found in previous work. In iKehrig et al.l yOOSj) a 
spatial separation between the WR star thereabouts and the Hen 
emission is found to be ~ 80 pc in the HII galaxy IIZwTO. More 
recently, IShirazi & BrinchmannI (120121) studying a sample of star- 
forming galaxies with nebular Hen emission from SDSS DR7 find 
that a large fract ion of the objects a t low metallicities do not show 
WR signatures. Kehrig et alj ( 120 llh also identify two Hen nebulae 
in ]y[33 which do not appear to be associated with any WR star. 

Different ionizing mechanisms other than WR stars (e.g. 
shocks; X-ray binaries) for the formatio n of the Hen line 
have been discussed i n th e literature (e.g., iGamett et al.l Il99ll ; 
IShirazi & BrinchmannI 12012 ). Shocks do not seem to be an im- 
portant excitation sour ce in Mrk 178 at large sinc e this galaxy is 
not detected in X-rays JStevens & Stricklandll998h , indicating that 
shock-ionization should not be responsible for the Hen emission. 
Our measured values of [Sii]/1/167 17,673 l/Ho- ratio (< 0.20) are 
lower than the ones obse rved in SNRs ([Sii]/1/16717,6731/Hq' ~ 0.5 
- 1.0; ISmithet alJI 19931] which points against shock ionization as 
well. The lack of X-ray detection also leads us to consider X-ray 
binaries not to be the main source of He+ ionization. 

Can we still associate the nebular Hen in Mrk 178 with its 
WR stars despite the spatial offset between the location of the 
WR stars and the Hen-emitting zones ? To answer this question 
we calculate the gravitational binding energy of a cloud with ra- 
dius of ~ 80-100 pc (~ the maximum observed distances between 
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Table 3. Results from the statistical analysis for the distributions of the physical conditions and chemical abundances in the ISM across the INTEGRAL FOV 
ofMrkl78. 



ISM property 



^weighted 



Statistical properties 



^weighted 



MGciiissicin 



CGai 



Sign.(%)'- 



log([On]/[Oin]) 


-0.56 


0.33 


— 


— 


<10 


n([Sii]) (cm-3) 


110 


95 


— 


— 


<10 


T([Oin]) (K) 


16770 


760 


17885 


2200 


26 


12+log(0/H) 
(direct method) 


7.72 


0.13 


7.69 


0.12 


80 


12+log(0/H) 
(N2 method) 


7.84 


0.16 


7.87 


0.16 


43 


log(N/0) 
(direct method) 


-1.40 


0.17 


-1.34 


0.16 


37 


log(N/0) 
(N202 method) 


-1.42 


0.19 


-1.38 


0.17 


75 



(a) en'or-weighted mean;(b) square root of the weighted sample vaiiance associated with the coiTesponding weighted mean; (c) mean of the Gaussian 
distribution; (d) standard deviation of the Gaussian distribution; (e) significance level of the null hypothesis in the Lilliefors test. We assume that the 
distribution of a certain property can be represented by a Gaussian fit if the null hypothesis is not rejected (i.e. Sig. > 10%) 
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Figure 10. log (iOm]A5007 /U/3) vs. log ([Nn]i6584/H(!'). Green circles 
show nebular Hen-emitting spaxels with no WR features, blue circles show 
the spaxels with blue WR bump (knot A in Fig.|5) and red circles the spax- 
els showing both the blue and red WR bumps (knots B and C in Fig.|5). The 
remaining spaxels are indicated with cross-symbols. Typical uncertainties 
on the line ratios are indicated by the cross on the top-left of the panel. The 
plotted measurements are corrected for reddening. The black solid curve 
is the theoretical maximum stai'burst model from Kewley et al. (2001), de- 
vised to isolate objects whose emission line ratios can be accounted for by 
the photoionization by massive stars (below and to the left of the curve) 
from those where some other source of ionization is required. 



the WR stars and nebular Hen across the FOV sampled here; see 
the map of Hen/14686 in Fig. [5) and assuming a hydrogen den- 
sity of ~ 100 cm"'. From this calculation we obtain that the en- 



ergy required to excavate a hole of this size is ~ 10^'"^^ erg. We 
checked that this number is of the order of the predicted mechan- 
ical energy injected into the ISM due to massive star winds tak- 
ing into account the number of WR stars derived in section 13.11 
These predictions are based on e v olutio nary synthesis models by 
iMolla. Garcia- Vargas. & BressanI (120091) for an HII region with 
Salpeter IMF and SMC-metallicity. Also, it is worthwhile remem- 
bering that WR stars with high temperatures (T* > 40 kK) and 
weak winds are e xpected to be capable of p roducing Hen ioniz- 
ing photons (e.g.. ICrowther & Hadfieldll2006l) . WC stars in knots 
B and C within Mrk 178 certainly appear to satisfy both criteria, 
although not exactly where nebular Heii4686 is seen, as discussed 
above. We should bear in mind that the precise geometry of the 
ionized ISM distribution in the close environment of these knots 
clearly plays a role. If WR stars from knot A were responsible for 
the nebular Hen emission in this knot, some of the WN stars would 
need to be early-types, which cannot be ruled out from our WR 
template fits analysis. Additionally, according to the spectral clas- 
sication scheme indicated in Fig. [TO] for most positions in Mrk 178 
our emission line ratios fall in the general locus of star-forming ob- 
jects suggesting that hot massive stars are the dominant ionizing 
source within our FOV. 



3.3 Aperture effects 

iBrinchmann. Kunth. & Durretl ( 120081) presented a study of WR 
galaxies based on spectra from SDSS DR6 and found Mrk 178 
to have the strongest WR features relative to Hj8 of any galaxy 
in the sample (see their Figs. 18 and 19) which is intriguing if 
one considers the low metallicity of Mrk 178. This curious be- 
haviour was possibly attributed to the fact that, due to the prox- 
imity of Mrk 178 (one of the closest systems in the sample of 
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1.5 2.0 

log EW(HP) 

Figure 11. EW(blue WR bump) vs EW(H/3). Asterisks show values ob- 
tained from SDSS DRV for metal-poor WR galaxies [7.0 <, 12-l-log(0/H) < 
8.1]; the red asterisk represents Mrk 178. The three blue circles, from the 
smallest to the biggest one, represent the 5. 7 and 10 arcsec-diameter aper- 
tures from our data centered at knot B. For the majoiity of the points the 
uncertainty of log[EW(Hj8)] is of the order of the symbol size. 



4 SUMMARY AND CONCLUSIONS 

We have presented here the first optical IFS study of the metal- 
poor WR galaxy Mrk 178. The proximity of Mrk 178 combined 
with the IFS technique allow us to locate and resolve star-forming 
knots hosting a few WR stars, and also to characterize the WR con- 
tent. In addition we are able to probe the spatial correlation between 
massive stars and the properties of the surrounding ISM. In the fol- 
lowing we list the main results derived from this work. 

• We defined three WR knots from which two are identified 
for the first time here. The WR knot spectra reveal the presence 
of nitrogen-type and carbon-type WR stars in Mrk 178. By 
comparing the observed spectra of the WR knots with SMC and 
LMC template WR stars we estimate a lower limit for the number 
of WR stars (> 20) in our Mrk 178 FOV that is already higher 
than that currently found in the literature (~2-3 WR stars from 
lOuseva. Izotov, & ThuanI 120001 . This is probably because IFS 
not only allows us to sample a larger area of the galaxy, but also 
offers a more power ful technique in sea rching for WR stars. As 
was demonstrated in iKehrig et al.l ( 120081 ) ■ using IFS one can find 
WR stars where they were not detected before because IFS, in 
comparison with long-slit spectroscopy, can increase the contrast 
of the WR bump emission against the galaxy continuum, thus 
minimizing the WR bump dilution. 



iBrinchmann. Kunth, & Durrell20o3) . the SDSS fiber could have by 
chance hit one single/few luminous WR stars while the Hfl flux is 
more extended. Actually looking at the Ha map in Fig. \5\ we see 
that the Balmer line emission is more extended and its peak is off- 
set from the WR emission which also indicates propagation of star 
formation; where the WR stars are now, the O stars are no longer. 
Here we take advantage of the IFS to investigate how the observa- 
tions of WR signatures are affected by aperture effects. 

In Fig.[n]we show EW(WR blue bump) vs EW(He) obtained 
from SDSS DR7 spectra of low metallicity [7.0 < 12-t-log(0/H) < 
8.1] WR galaxies and we find the s ame discrepancy reported by 
IBrinchmann. Kunth, & Durreu (|2008J), i.e. the most deviant point 
(red asterisk) in Fig. [TT] belongs to the galaxy Mrk 178. Using 
our IFU data we construct one-dimensional spectra by combin- 
ing spaxels within circular apertures of increasing diameter (the 
center of all apertures coincides with our WR Knot B which is 
approximately the center of the Mrk 178 SDSS fiber). The mea- 
surements derived from these spectra are blue circles in Fig. 1111 
which show that Mrk 178 is brought closer to the bulk of metal- 
poor systems when the aperture size increases confirming that the 
offset between Mrk 178 and the rest of the objects as viewed by 
SDSS is not real and caused by aperture effect instead. We also 
checked that for apertures with diameter larger than ~ 10 arc- 
sec we no longer detect the WR bump. This likely happens be- 
cause the underlying stellar absorption masks the WR stellar sig- 
natures more easily if the apertures are too large. In metal- poor 
systems there is an enhancement of this effect since WR stars are 
expected to have weaker lines than in more metal rich environm ents 
dConti, Garmanv. & Massevlll989l : ICrowther & Hadfieldll2006l) . In 
addition, after correcting the EW(WR blue bump) by the physical 
area of the SDSS aperture for each object in Fig. [TT] we verified 
that the measurements of EW(WR blue bump) for all galaxies (not 
only Mrk 178) suffer from so me degree of aperture effects (see also 
IPerez-Montero & Diazll2007l) . All this can explain why sometimes 
WR feat ures are not detected when they are expected to be present 
(see also ljames, Tsamis, & BarlowiuOlOl) . 



• Our statistical analysis suggests that spatial variations in the 
gaseous electron temperature exist and that the scatter in Tg can 
be larger than that in 0/H within the observed FOV. Thus, caution 
should be exercised when analysing integrated spectra of giant HII 
regions/HII galaxies which do not necessarily represent the "local" 
ISM properties around massive star clusters. 

• The nebular chemical abundance in Mrk 178 is homogeneous 
over spatial scales of hundreds of parsecs. This result is in 
agreement with previous work on ionized gas in HII ga laxies and 
favours the scheme, discussed by iTenorio-Tagld (119961) , in which 
the most recent products of massive star nucleosynthesis remain 
unobservable optically on timescales longer than the lifetime of 
HII regions. The representative metallicity of Mrk 178 derived 
here is I2-l-log(0/H)=7.72 ± 0.01 (error-weighted mean value of 
0/H and its corresponding statistical error). 

• To probe the presence of small-scale (~ 20 pc, our resolution 
element size) localized chemical variations, we performed a close 
inspection of the chemical abundances for the WR knots from 
which we find a possible localized N and He enrichment, spatially 
correlated with WR knot C. 

• The nebular Hen emission appears extended and is likely 
related to ionization from hot stellar continua. We interpret 
the spatial offset between WR stars and nebular Heii-emitting 
regions as an effect of the mechanical energy injected by WR 
star w inds This agrees with results by Shirazi & Brin chmannI 
( 120 12h and iKehrig et alj ( 120081) . Caution should be exercised 
when the precise geometry of the ionized ISM is not considered. 
Shock-ionization and X-ray binaries are unlikely to be signif- 
icant ionizing mechanisms since Mrk 178 is not detected in X-rays. 



• IBrinchmann, Kunth, & Durretl ( 1200 8h report an unusual (too 
high) L(WR bump)/L(I^) ratio for Mrk 178 derived from its SDSS 
spectrum. The authors question this assuming that aperture effect 
can be an explanation. This is confirmed in this work with IFS. We 
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also demonstrate that using too large an aperture, the chance of de- 
tecting WR features decreases. This result indicates that WR galaxy 
samples/catalogues constructed on single fiber/long-slit spectrum 
basis may be biased in the sense that WR signatures can escape de- 
tection depending on the distance of the object and on the aperture 
size. 

This paper shows the need to identify good targets, such as 
Mrk 178, and the power of IPS for investigating issues related with 
aperture effects. In particular, we show the effect of distance and 
metallicity on the search for WR features, and the extent to which 
physical variations in the ISM properties can be detected. In addi- 
tion, we believe that this work makes way for future research of low 
metallicity, close-by objects aiming to survey and characterize the 
WR population and to identify Hen nebulae through high-spatial 
resolution imaging and IPS. 
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